par(mfrow=c(2, 2))
data1_16 <- read.csv("data/时间序列分析——基于R（第2版）案例数据/csv/A1_16.csv")
ts_x <- ts(data1_16$x, frequency = 4, start=1981)

# 时序图
plot(ts_x, type="o", pch=5, col='#39CBB4')

# 进行 Holt-Winters三参数指数平滑（乘法模型）
x_fit2 <- HoltWinters(ts_x)
x_fit2
# 进行 Holt-Winters三参数指数平滑，进行12期预测
# install.packages("forecast")
library(forecast)
x_fore <- forecast(x_fit2, h=12)
x_fore
# 预测效果图
plot(x_fore, lty=2)
lines(x_fore$fitted, col=4)
